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Abstract 

Recently it has been experimentally demonstrated that certain glasses display an unexpected 
magnetic field dependence of the dielectric constant. In particular, the echo technique experiments 
have shown that the echo amplitude depends on the magnetic field. The analysis of these experi- 
ments results in the conclusion that the effect seems to be related to the nuclear degrees of freedom 
of tunneling systems. The interactions of a nuclear quadrupole electrical moment with the crys- 
tal field and of a nuclear magnetic moment with magnetic field transform the two-level tunneling 
systems inherent in amorphous dielectrics into many- level tunneling systems. The fact that these 
features show up at temperatures T < WOmK , where the properties of amorphous materials are 
governed by the long-range R~ 3 interaction between tunneling systems, suggests that this inter- 
action is responsible for the magnetic field dependent relaxation. We have developed a theory of 
many-body relaxation in an ensemble of interacting many-level tunneling systems and show that 
the relaxation rate is controlled by the magnetic field. The results obtained correlate with the 
available experimental data. Our approach strongly supports the idea that the nuclear quadrupole 
interaction is just the key for understanding the unusual behavior of glasses in a magnetic field. 

PACS numbers: 61.43.Fs, 64.90. +b, 77.22. Ch 
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I. INTRODUCTION 



In 1998 Strehlow et al.— observed a pronounced but unforeseen magnetic field dependence 
of the dielectric constant of certain multicomponent glasses at ultralow temperatures T, i.e., 
below lOmK, for fields as small as 10/iT. The effect was found to be by several orders 
of magnitude larger than expected for an insulator in the absence of magnetic impurities. 
This result was especially astonishing since earlier experiments did not indicate noticeable 
magnetic field effects 2 *. Yet, careful measurements^ in magnetic fields ranging up to 25T 
and temperatures below lOOmK revealed that the magnetic field causes drastic changes in 
the dielectric response. These experiments have caused further investigations 4 -"*^. Follow- 
ing the original discovery, a number of different properties of the glasses were investigated 
in magnetic fields at these temperature o 7 ' 8 ' 9 ' 10 ' 11 , e.g., by using the dielectric dipole echo 
technique 1 ^. 

Many low-temperature properties of glasses have been successfully described in the past 
by a standard tunneling mode l 13 ' 14 ' 15 . To a good approximation a Tunneling System (TS) can 
be treated as a particle moving in a double- well potential (DWP). Because of the randomness 
of the glassy structure, the energy difference between the two wells as well as the tunneling 
matrix element have a broad distribution. This distribution is practically universal for all 
known dielectric glasses and results in an agreement between theory and experiment above 
100m.fr. Yet, below this temperature it is necessary to extend the model of isolated TS's by 
taking into account the long-range interaction between them in order to interpret numerous 
experiment a 16 ' 17 ' 18 ' 19 ' 20 ' 21 . The concept of resonant pairs (RP) plays an important role here. 
It captures the important physics as long as the interaction between TS is not so strong that 
the tunneling picture looses its meaning. 

After the discovery of the anomalous glass behavior in a magnetic field, several extensions 
of the standard TM have been suggeste d 22 ' 23 ' 24 . The dielectric properties of glasses at low 
temperatures are known to be due to the TS. It is reasonable to suppose that this is also the 
case for glasses in magnetic fields. Therefore, the principal question is how the magnetic field 
is acting on the TS's. In our opinion, the model of Wiirger, Fleischmann and Enss 2 ^ deserves 
special mentioning and attention. In the framework of the model a direct coupling between 
the nuclear spin of a TS and the magnetic field takes place, though initially it was alleged that 
this possibility should be ruled out 4 -. However, the echo experiments convincingly evidence 
the influence of the nuclear moments on tunnelin g 7 ' 8 ' 9 ' 10 . Recently Nagel et al.— investigated 
isotope effects in polarization echo experiments in glasses. They observed magnetic field 
effects on amorphous glycerol when hydrogen, which has no electric quadrupole moment, 
was substituted by deuterium, that possesses nonzero quadrupole moment. Thereupon, one 
can conclude that the quadrupole electric moment of the TS's is the key feature responsible 
for the effect. 

A generalization of the standard tunneling model can be done as follows 2 ^. Consider 
a tunneling particle with a nuclear spin. The energy levels of the system are degenerate 
with respect to the nuclear spin projection. This degeneracy is split if the particle has a 
quadrupole electric moment Q which interacts with an electric field gradient (EFG), i.e., 
qij = d 2 V/d 2 XiXj, where V(r) is the crystal-field potential. For the sake of simplicity, 
Wiirger et. al.— supposed that the EFG possesses axial symmetry. In this case, the energy 
of the tunneling particle depends on the orientation of the quadrupole quantization axis. In 
general, in glasses the axes and differ in the two wells of the potential (see FigC]). For 
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Figure 1: Tunneling system with asymmetry energy A, and the corresponding quadrupole quanti- 



this reason, the quadrupole energy changes when the particle tunnels through the barrier. 
The magnetic field influences the energy spectrum since the particle acquires a Zeeman 
energy depending on the nuclear spin projection. Within that generalized tunneling model 
the previously mentioned echo experiment could be described^. 

The unusual dielectric response of glasses to an applied magnetic field occurs in the 
temperature region T < 50mA" where the standard tunneling model of non-interacting TLS 
cannot explain many experiments 17 . Just within this temperature region RP are responsible 
for the relaxation in glasses. For this reason it is interesting to investigate whether or not 
the application of the magnetic field influences the dynamics of resonant pairs. 

The principal goal of the paper is to show how the quadrupole interaction, if any, shows 
up in different characteristics of glasses. For this purpose we discuss first the explanation of 
the anomalous magnetic-field induced behavior of dielectric properties of glasses at ultra-low 
temperature based on the model^i Then, we clarify how these features are reflected in the 
relaxation phenomena governed by the interaction between tunneling systems. We shall 
extend in the paper the concept of RP so that it includes the quadrupole interaction as well 
as the interaction with the magnetic field. We will also investigate how theses interactions 
govern the relaxation rate of tunneling systems. A comparison of some of the results obtained 
with experiments will be made. 

II. GENERALIZATION OF THE STANDARD TUNNELING MODEL 

Consider a double-well potential (DWP) characterized by the asymmetry energy A. At 
low enough temperatures only two energy levels corresponding to the ground state in each 
well of the DWP are significant. These levels are connected by the tunneling amplitude A . 
According to 13 i 14 i 15 , the parameters A, A obey the universal distribution 



where P as a constant. An isolated tunneling particle is usually described by the standard 
two- level pseudospin 1/2 Hamiltonian 



zation axes in the left well ul and in the right well ur (see RefJ&) 




(1) 



h 



— a cr . 

2 2 



(2) 



3 



L 

a) b) 
Figure 2: Energy levels: (a) in absence and (b) in the presence of an applied magnetic field 



Suppose that the tunneling particle possesses nonzero spin I 2 = 1(1 + 1). Below we 
will identify it with the nuclear spin. Then the states of the particle are characterized by 
the sign of the pseudospin projection and by the particle spin projection onto a proper 
quantization axis n = —I, Thus, the dimension of the Hilbert space for the tunneling 
particle is 2 (21 + 1). Transitions of the particle between the wells of the DWP occur with 
conservation of the spin projection. Introducing the spin projection operator \n) (n\, one 
can define generalized pseudospin operators o z [n] and a x [n] so that a 1 [n] = a % ® \n) (n|. 
As long as the spin of the particle does not interact with the environment, the energy levels 
are degenerate with respect to the spin projection and the tunneling Hamiltonian reads 

fc = -W£a>] + ^*[„]Y (3) 

In this case, the spin degrees of freedom do not exert an influence on the tunneling properties 
of the particle, no matter which Hamiltonian, either (j2J) or (jSJ), is used to describe the particle 
motion. 

Suppose that an uniform magnetic field B directed along the z— axes is applied. Then, 
the tunneling particle gains extra Zeeman energy dependent on the spin projection I z 

E int = g(3BI z , (4) 

and the degeneracy of the energy levels is lifted. However, this splitting is irrelevant. Indeed, 
in both wells of the DWP the magnetic field has the same magnitude. For this reason, the 
Zeeman contribution depends only on the spin projection and does not depend on the 
pseudospin projection. For the case 1 = 1, the energy structure of the tunneling particle 
before and after application of the magnetic field is presented in Fig. (j2J). The states with 
fixed spin projection are the eigenstates. One should pay attention to the fact that the 
tunneling between the two sites L and R can happen only between eigenstates that have 
equal spin projection (like in the absence of a magnetic field). So, the magnetic field does 
not influence the overlap integral between the wave function of the left and the right well. 
This means that the application of a magnetic field alone does not influence the properties 
of the TS under consideration. 

Let us now consider the case when the spin of the tunneling particle is I > 1. Then 
the tunneling particle possesses a quadrupole electric moment. It interacts with the crystal 
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field which is characterized by the electric field gradient (EFG) q^. The Hamiltonian of the 
particle interacting with the crystal field reads^ 

Hq = - 2I( lf_ 1) [9iiA 2 + Haft + fe/l] (5) 

where Hq is written in the basis e±, e2, in which the tensor has a diagonal form. The 
EFG satisfies Laplace's equation q u + g 2 2 + <?33 = 0. Introducing the asymmetry parameter 

x = — ; > ( 6 ) 

Qu 

one can rewrite Eq. (JSJ in the form 

ff = -6(3i? + x(i|-i|)-I a ). (7) 

Here the parameter 6 = 4/^"-^ designates the quadrupole interaction constant. We assume 
that the ei, e^, e3 axes are chosen so that #33 < ^22 < since then < x < 1. If x = 0, the 
EFG possesses axial symmetry. In this case the quadrupole energy is completely defined by 
the spin projection Ii and the quadrupole quantization axis is directed along e x . 

Assume that in the left well the basis is e\^e^e^ while in the right well it is e'^e^e^. 
In general, these basises are different. Let us introduce a basis e x , e y , e z common for both 
wells. We assume that e\ and e[ lie in the e x ,e y plane. Suppose that e\ coincides with e x 
while e[ forms an angle 6 with the e x axis. Then, for the case 1 = 1 considered here, one 
can represent Hq in the right well in the basis of the eigenfunctions of the operator I z , 
I — X> , 1 0> , I X> , as follows* 

/ -J(l + x) |(l-f)e-" 
H Q (0) = b\ (x+1) j. (8) 

V 1 C 1 - f ) e+ie -K 1 + ^) 

Hq (6) has the following eigenstates and eigenvalues 




(9) 



For the left well the eigenvalues are the same as for the right well while the eigenvectors 
are \a°) , {a 1 ) , |a 2 ), resulting from Eq.Q for 6 = 0. It follows directly from (JSJ) that the 
spin state \a°) = |/3°) = |0) is orthogonal to all the others, but because of the fact that the 
quadrupole quantization axes differ in the two wells, (a l \ ft) 7^ 0, i ^ j. Thus, the state 
with spinor |0) in the left well can couple only to the state with spinor |0) in the right well, 
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Figure 3: Energy levels structure induced by the quadrupole splitting and possible tunneling transi- 
tions between them. In zero magnetic field the distance between the levels with \a°) and [ct 1 ) , and 
the levels with \a°) and \(3i) is defined by the parameter k and vanishes as this parameter goes 
to zero. 



while the remaining states in the left and the right well couple to each other. The energy 
level structure and possible transitions between the left and the right well states are shown 
in Fig. ©. 

Let us investigate how the application of the external magnetic field changes the energy 
spectrum of the tunneling particle. The magnetic field shows up in two ways. First, it 
produces a Zeeman splitting. Then for different TS's the planes formed by the vectors e%, e[ 
are randomly oriented. For this reason the magnetic field is in general not orthogonal to this 
plane. Therefore it mixes the states with spinor |0) with states with spinors \a l ) and \ft) 
(i = 1, 2) of the two wells. In a new basis \a l ) and \ft) which is dependent on the direction 
of the magnetic field all (a l \ ft) ^ 0, i,j = 0,1,2. Thus, additional transitions between 
the two wells become possible (see Fig. (jlj)). 



Figure 4: Shown in the figure by thick lines are the extra transitions created by the application of 
the magnetic field. The thin lines correspond to the previous transitions (see Fig. . 
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It should be stressed that with increasing magnetic field, the Zeeman splitting will even- 
tually exceed the quadrupolar one. Therefore, we discuss first the case when the TS's are 
affected only by the magnetic field (see Fig. (J2J). 

To simplify the further analysis, suppose that the magnetic field is orthogonal to the 
plane ei,e[ and that the quadrupole splitting constant b is the same in both wells. Then 
the Hamiltonian of the TS in the presence of the quadrupole and Zeeman splitting reads 



H L A -I 
A -I H R 



(10) 



where I is the rank 3 unit matrix and 



A - + |(l-f) 



Hl= ( 2 |(1 -f) -# + 6(x+l) 2 3 ] (11) 

-f + m-|(l + x) 



| - m - | (1 + x) | (1 - f ) e- w 

H R = | f + 6(x+l) | . (12) 



\(\-%)e id ~ #+m-!(.i .-rx) 



ILL PERMITTIVITY OF A TUNNELING SYSTEM WITH QUADRUPOLE AND 
ZEEMAN SPLITTING 

Let us investigate how the changes of the energy spectrum induced by the quadrupole 
and Zeeman splitting described by Hamiltonian (jl())l influence the properties of a TS. First 
we examine the influence of these changes on the dielectric permittivity. 

Consider a particle that can occupy {21 + 1) levels in a DWP. We introduce the dipole 
moment operator in the form 

p=j[o I ) ^ 

where I is the unit matrix of rank (21 + 1) and po defines the value of the dipole moment. 
(For 1 = this expression transforms into the well known dipole moment operator for a 
two- level-system) . The interaction of the TS with the external electrical field F reads 

V = -Fp . (14) 

In second-order perturbation theory the correction to the energy of a TS induced by this 
interaction is given by 

Se=F > z -^Lmmi (15) 

b=f=a 

where T is the temperature, Z = ^2e~^r~ is the partition function and E a denotes the 

a 

eigenvalues of the TS in the presence of a magnetic field. 
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On the other hand, as the electric field changes from zero to a certain final value F, the 
energy increases by X-^ 2 /2, where \ is a permittivity of the system. Comparing this value 
to be one obtains 



X 



2Z 



"'E 



e t 



* |(a|p|6)| 5 



■'E 



r — e 



|(a|p|6)| 



-Eh — 



(16) 



Let us apply this relation for the TS described by Hamiltonian (jl(Jj) . 

The permittivity of a tunneling system is a function of the parameters A, A , b, 8, m. 
(For a preliminary analysis in this part of the paper we confine ourselves to the case of 
an axial symmetric EFG by setting the parameter x = 0). Keeping in mind a possible 
fit of the available experimental data, we choose the following values for the parameters 
A max = Aomax = 10 K, b = lOmK. The temperature varies between lOmK and lOOmK. The 
Zeeman energy of the nuclear spin ranges from zero to 30mK. With the Lande factor g = 3, 
this corresponds to a maximal magnetic field of about 15T. 

The permittivity has been estimated numerically. First, we calculated the eigenvalues and 
eigenfunctions for the matrix (jl(Jj) with values of the TS parameters in the above mentioned 
range. Then, we estimated the permittivity \ (A, A , b, 8, m, T) by using Eq. (|16Jl . Finally, 
the result was averaged over the parameters A, A , 8 by using Eq.Q and assuming a uniform 
distribution of the angles formed by the quadrupole quantization axes in the wells. As an 
example, we present here the result of our calculations for the temperature T = 40mk 




Magnetic Field B [T] 



Figure 5: Magnetic field dependence of the permittivity e' for typical quadrupole splitting b = WmK 
atT = AOmK 



The permittivity exhibits a pronounced peak at an energy approximately corresponding 
to the magnetic field B w 15T. This result corresponds well to Ref.— where a pronounced peak 
in the permittivity was observed for a similar value of the magnetic field and a temperature 
of T = UmK. 

As mentioned before we have assumed x = in calculating the permittivity shown in 
Fig. El When we introduce x a second energy scale appears when x is small. This scale 
is relevant for weak magnetic fields. It might explain another peak observed in Ref.— in the 
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low field regime. Also a second atom of the tunneling entity with a different quadrupole 
moment may result in a second energy scale. We have modelled a low energy scale by simply 
redoing the calculations for b = 0.3mK and plotting the results for the permittivity in a 
form which can directly compared with the experimental findings of Ref A Figure |H1 shows 
an agreement with experiments that the permittivity at fixed value of magnetic field is the 
higher the lower the temperature is. 
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Figure 6: Magnetic field variation Ae of the permittivity. 



IV. RELATION BETWEEN THE PERMITTIVITY AND MANY-BODY RELAX- 
ATION 

As mentioned in the introduction, the temperature region T < lOOmK where certain 
glasses show an unusual response to an applied magnetic field coincides with the region 
where the relaxation of tunneling system is due to the R~ 3 interaction between TLS rather 
then single phonon processes. In the previous section we have seen that a magnetic field 
affects the permittivity provided the quadrupole effect is taken into account. Therefore, it is 
of interest to investigate whether a magnetic field influences the relaxation induced by this 
interaction. 

The relaxation induced by the long range R~ 3 interaction is strongly connected with the 
concentration of resonant pairs (RP). Resonant pairs of TS's are the main concept for this 
relaxation mechanism. First we recall briefly the main idea of this approach. 

Consider a pair of two-level-systems. In general, it possesses four different configurations. 
Only two of them (configurations A and B ) shown in Fig. (J7J) are important. The special 
feature of these configuration is that one TLS is in the ground state while the other is in 
the excited state. 

Let V\2 = U/R 3 denote the weak interaction between these TLS's where U is the inter- 
action constant. Such a pair is in resonance when 

\E X -E 2 \< V 12 . (17) 
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Figure 7: The states A ) and B ) belong to the flip - flop configuration and form a resonant pair. 
The states C) and D) are energy-wise forbidden 



Here E\ = \J A^ + Af, Aoi, Ai,i = 1,2 are the tunneling parameters of two TLS constituting 
a pairs. We are interested in the case E\ « E 2 ~ T and Vi 2 <^ T which means a weak 
interaction. Due to the constraint (|17|) . the states (A) and (B) are coupled with each other 
and are well separated from the states (C) and (D) by an energy gap of the order of T. For 
this reason, they can be excluded from the consideration. The two states of the TLS pair (A) 
and (B) are referred to as a flip-flop configuration and the transition between them is called 
a flip-flop one. When such two TLS's form a resonant pair (RP) they can be considered as 
a new type of the two - level system with the asymmetry | E\ — E% \. Resonant pairs are 
responsible for the many-body relaxation induced by the R~ 3 interaction. The transition 
amplitude between the levels (A) and (B) is described by the following expressio) 



,16,18 



17 
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Vi 2 a^a. 



x AoiA()2 

E\Ei 



(18) 



Next we write down the expression for the concentration of pairs with tunneling parameters 

A A 16.18 . 



p(2) ( A 



pi Ao P ) 



P 2 I 



Aoi 



d A 02 dA 2 
A 02 



l + e~ 



1 , ^ 
1 + e t 



J d*R5 (A 0p - ^ §) 5 (A p - (E, - E 2 )) 



(19) 



The integration over d 3 R gives 



d s R5 



A 0p 



AqiA 02 U 
E X E 2 R 3 



A01A02 U 



E\E 2 



A 2 ' 



(20) 



T and, therefore, one can omit A p in 



For resonant pairs one has A p < Ao p -C E\ « E 2 

the argument of the 5— function in ()19j) . The concentration of resonant pairs Pr 2 ^ (Aq p ] 



is 
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obtained by integration of (A p , An p ) over the interval < A p < An p 



P?) (Aop) , P ^j^. s (^) , {El _ *, (21) 

On the other hand, the resonance permittivity of a TLS is given by the well known 

14 \ 



expression [see e.g.- 



B 



dA 1 /AoY . E 



o JSo An £? V E J 2T 



X = P dA ^ tanh- (22) 



Let us calculate the derivative 



* ^ B dA r^ifv, (23) 



dlnT 2T7 4 A c/i 2 f 

Comparing expressions (}2~T]) and (J23j) . we can conclude that P r (A 0p ) and the square of 
d ^ T are approximately proportional to each other. Therefore we expect similar features 
in the behavior of the permittivity and in the rate of the interaction induced relaxation 
which is proportional to the RP concentration. This means that the multilevel systems 
can again be considered as effective two level systems. When the Zeeman splitting exceeds 
the quadrupole one, the multilevel tunneling systems behaves exactly as two-level-systems. 
In Figure |S] we present in the framework of the model described above the results of the 
calculation of the parameter ( 9 ®* T ) for different temperatures as a function of the applied 
magnetic field. One notices that for magnetic fields causing a Zeeman splitting smaller or of 
the order of the quadrupole splitting, i.e., b < lOmK it holds that the larger the temperature 
the larger is dx/dhaT. This can be understood as follows: The relaxation rate induced by 
the interaction between TS's is proportional to the concentration of RP, which in turn is 
proportional to the temperature. This conclusion agrees with the sequence of different curves 
in Fig. [HI However, some curves show a monotonous magnetic field dependence, while others 
display a non-monotone behavior. This feature is explained in the following section. 



V. RELAXATION OF THE INTERACTING MANY-LEVEL TUNNELING SYS- 
TEMS 

The main goal of the current section is a qualitative treatment of the results obtained in 
Sec. IV. In particular, we want to explain the temperature dependencies shown by Fig. |H1 

In our previous paper n 16 i 18 i 19 i 20 i 21 we have demonstrated that the many-body R~ 3 inter- 
action between TLS's results in a new relaxation mechanism responsible for the low tem- 
perature relaxation. On the other hand, it was shown above that if the tunneling particle 
possesses a nuclear spin / > 1, the energy spectrum of a tunneling system consists of several 
lines, as distinct from one line in the case of a TLS. Below we investigate how these changes 
in the energy spectrum of the tunneling system caused by the interaction of the tunneling 
particle with the electric field gradient and the magnetic field influence the relaxation rate 
produced by the R~ 3 interaction. 

The model under investigation assumes that the particle can occupy n = 21 + 1 levels 
each in the left and in the right well of a TS. Let us generalize the concept of resonant pairs 
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Figure 8: The square of the logarithmic derivative of the permittivity as a function of the magnetic 
field. 
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Figure 9: Two tunneling systems:TS 1) is in one of its ground state multiplet. TS 2) is in the 
excited state multiplet. 

The occupied states are marked by solid line. The arrow lines show which levels are occupied after 
flip-flop process has happened 



12 



to the case n > 1. Consider the tunneling system (1) in Fig. El There we are dealing with 
a ground-state multiplet the states of which are denoted by their spinor parts a\ and an 
excited state multiplet denoted by (3[. The corresponding energies are 5 a ^ and 5^. They 

are derived from the energies and E\ = a/ A 2 + A 2 of the TS in the absence of the 
quadrupolar and Zeeman interaction. Therefore the various transition energies between the 
two multiplets are 



5 



hJ 



n 



(24) 



Like in the case of a TLS only those with E[ J 
processes. Below we will assume that also 



T contribute significantly to relaxation 



5 a i , 5 £ 



< T. 



(25) 



holds 31 . Therefore only E\ defines the Gibbs distribution and the states within a multiplet 
are occupied with nearly equal probability. A similar analysis holds for the TS(2) of Fig. 
El Thus, two multilevel TS's give raise to flip-flop processes a\, /3 2 =>■ (3± , a l 2 when the 
corresponding energy differences between the two configurations differ in energy by less than 
a characteristic interaction matrix element. The latter and the various different processes 
are discussed in the following. 

Let V [a\, f3 2 ; P%i a %\ denote the transition amplitude of such a flip-flop transition. Then 
the condition for a resonant pair is 



kl\ 



E 1 -E 2 + [5 



\E?-E 



When all the energy splittings S 
given by 16 i 18 i 19 i 2 °i 21 



5 P{ 



8/3* 



< 



V[a\,(3i;^,a 2 ] . (26) 



0, i.e., for TLS's the resonating pair concentration is 



N* = (PT) (PU) 



(27) 



Let us investigate how the concentration of resonating pairs is modified for multilevel 
systems. For simplicity we set all V [a\, P 2 ', Pi, a l\ — V- We fix the energy E\ and also 
the splitting energies 5™. Then the energy range of the parameter E 2 is defined by Eq. 
fl2TI|) and equals V. For a given E\ there are n A different state configurations a\, a 2 , f3° 2 

possible, i.e., the energy difference — 8^ — {& a \ ~ takes n 4 different values. So 

there are also n A different energy ranges for E 2 . Different values of E 2 relate to different 
tunneling systems. However, one must take into account that the probability of finding an 
initial configuration a\, (3{ is 1/n 2 . Hence the total probability of forming resonating pairs 
would seem to increase like n A /n 2 = n 2 as compared with TLS's. This would hold true 
if the transition amplitude V [a\, {3 2 , P{, a 2 ] would coincide with the previously discussed 
amplitude Vi 2 of TLS's. But this is not so as is easily seen as follows. The transition between 
two TS's (1) and (2) can be described similarly as in Eq. (|18j) by 



(28) 



V [a\,Pi;Pi,a 2 ] = V X2 {P{\ a\) {a 2 \ ft) 



(29) 
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As before the notation \a\) , \ etc. refers to the spinor part only of the wavefunction 
of the particle. Here | V12 1 = A J^° 2 £//-Ri 2 an d -R12 is the distance between the TS. Let 
us estimate the transition amplitude V [a\, /3f, a\\ in terms of V\ 2 - By making use of 

■ 1 2 

Eq.(j2HI) and the completeness condition in the form of a \)\ = 1> we obtain 



|^H^2;^,^]r = i^ 2 | 2 ^K^|ai>rEK^I« 2 >| 



= \V 12 \ Z (30) 

For further estimation, we make the simplest assumption that all the matrix elements (f2§|) 
entering Eq. (J30)) are equal. The total number of terms in the left-hand sum is n 2 . So 

\V[a\,Pi;Pi,ai]\ = ^- (31) 

Thus, the tunneling amplitude V [a\, (3{; a 2 , (3 2 ] entering Eq. (ffilj) is n times smaller as 
compared with the case of pairs formed by TLS. Therefore, the factor n 2 found above for 
increasing the probability of formation of resonating pairs is further reduced by a factor n. 
So, the total probability to form a resonating pair increases for multilevel systems £ = n 
times as compared with TLS's. 

A similar analysis can by carried out when some of the energy levels remain degenerate 
or when the transitions between some of them are forbidden. This case is investigated by 
using as an example the spectrum of TS with a nuclear spin 1 = 1 when only quadrupolar 
effects are taken into account (see Fig. EJ). In order to find the matrix element of the flip-flop 
transition let us again make use of Eq. ()29|) . Here we must differentiate between three cases. 

In the first case, the transition in both TS's of the pair occurs between levels with I z = 0. 
The corresponding states are denoted by their spinor parts a?, f3® and a®, ft®- By taking into 
account that the spinor parts remain unchanged it follows that 

V [a?, a°] = V 12 <a?| $) a° 2 ) = V 12 . (32) 

In the second case, the transition between I z = states takes place only in one TS of the 
pair. In that case one finds 

V[a 1 ,Pi-,(3° 1 ,aj i ]=V 12 (a 1 \p 1 )(Pi\aj 2 ) . (33) 

The sum of the square of the matrix elements is 

\V K&fooQ f = \V l2 \ 2 J2 M «a>f = l^l 2 • (34) 

■ 1 2 

Hereby the condition liPl] a 2)\ = 1 nas been used. Since the total number of terms in 



the sum (|3*3j) is two, one can estimate 

Rl- R9 I = . 

V2 



7[a?,^;^,c4]| = C^l . (35) 
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In the third case, when the transitions are between levels with nonzero spin projection in 
both TS's, one finds 

|^K,^;^,^]| 2 = i^ 2 | 2 ^|<^>orEK^I«2>r 



« 



= \V 12 \ 2 . (36) 
The total number of terms in the left-hand sum in Eq. (}3l)j) is equal to 2 2 and, therefore, 

\V[a\,Pi;Pi,ai]\ = ^ . (37) 

Using Eqs. (|2l)j) . (J32J) . . (}3T|) one estimates that the probability to form a resonating pair 
increases by a factor of £ = 1 + ~ 1.6. 

Thus, for the case of a nuclear spin 1 = 1 the quadrupolar interactions result in an 
increase of the probability of finding resonating pairs by a factor £ « 1.6 as compared 
with TLS's. When in an applied magnetic field the Zeeman splitting is of the order of the 
quadrupolar interaction, the probability of forming a resonating pair increases by a factor 
of £ = 3. Finally, when the Zeeman energy exceeds the quadrupolar splitting, the factor is 
£ = 1. 

The results obtained in this section are based on the fact that a tunneling system is a 
multilevel one. In other words, the energy levels of the TS should be well resolved. Yet, in 
an ensemble of interacting TS's the energy levels fluctuate due to spectral diffusion^. When 
the scale of spectral diffusion exceeds the quadrupole splitting the transition occurring from 
the different levels of the TS's can not be considered as statistically independent. The scale 
of spectral diffusion is about 7T, 7 = PU—. So our approach is valid when the quadrupole 
splitting b is 

b > 7 T . (38) 
Due to a similar reason the Zeeman splitting should obey 

g f iB> 1 T . (39) 

This condition establishes the minimal value of the B in order that magnetic field effects 
show up. 

The results obtained in this section are based on the energy level classification described 
in Sec. II. This requires that the tunneling amplitude A fulfills the relation A -C A, b. For 
that reason, in the above analysis the factor ^^° 2 entering Eq. (|18|) is a small parameter 
formed by strongly asymmetric tunneling systems. 



VI. DISCUSSION AND CONCLUSIONS 

We have shown that the strong magnetic field dependence of the electric susceptibility 
in ultracold glasses can be understood by taking into account the interactions of tunneling 
systems in the presence of nuclear quadrupolar moment. The essential point is the following: 
Even by a small applied magnetic field the number of different energy levels of a tunneling 
system is increased. This in turn, modifies the concentration of resonant tunneling pairs 
and leads this way to observable effects. 
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Two kinds of experiments can be related to the present investigation. The first kind deals 
with the experimental determination of the real part of the response of the system to an 
applied magnetic field. We find for the field dependence of the static electric susceptibility 
the right order of magnitude of the effect. That a too small effect was found in^ 2 - may result 
from the use of perturbation theory. Another unexplained feature, namely a plateau in the 
temperature dependence of the electric susceptibility sometimes called dielectric saturation 
is most probably due to a dependence of the tunneling matrix element A on the spinors of 
the right and left well. 

The second group of experiments is related to measurements of the imaginary part of 
the response. In this case information on the relaxation rate of elementary excitations can 
be obtained. For example, measuring the echo amplitude allows to determine the behavior 
of the TS coherence time 72 as function of field. The echo amplitude has been found to 
show a non-monotonic dependence on the magnetic fiel d 5 ! 6 ! 7 ! 8 ! 9 ' 10 ' 11 . It is defined by the 
transverse relaxation rate r^ 1 . In Sec. 5 we have introduced the parameter £ depending on 
the magnetic field which controls the concentration of RP. If the relaxation in the system 
is due to the resonant pairs, the transverse relaxation rate is directly proportional to the 
concentration of RP. Without the quadrupolar effect, the relaxation rate r^ 1 = 7 2 T— . If 
the quadrupolar effects are taken into account, the concentration of RP increases by the 
factor £. Therefore, the relaxation rate becomes 

T^ntfT . (40) 

The non-monotonic behavior of the parameter £ on the magnetic field correlates with the 
behavior of the echo amplitude found in the experiments mentioned above. Note that in 
high magnetic fields £ reduces to £ = 1, because in that limit the energy levels are equally 
spaced as in Fig (J2J). 

The standard scheme of interpreting the echo experiments is based on the assumption 
that TS's are two-level ones. In that case they can be described in terms of r^ 1 . This 
is not the case when the quadrupolar effects are taken into account. Therefore it is no 
surprise that the calculated r^ 1 can only qualitatively but not quantitatively describe the 
echo experiments. Instead the approach i n 24 i 34 i 35 based on a multilevel description of TS's 
instead of the standard Bloch equations seems very reasonable. 

For this reason, another kind of low-temperature experiments should be made to clarify 
the role if the R~ 3 interaction between the tunneling system. For example, measurements 
of the magnetic field dependence of dielectric loss or internal friction would allow to extract 
information on the lifetime T\ of the excitations. We want to show that in this case the 
factor £ should show up in an even more important way. 

The spectral diffusion reaches its maximal value 7T at time T\. Therefore, the rate of the 
spectral diffusion is jT/ti. On the other hand the spectral diffusion rate relates to r ~ 2 18,27 ^ 
Then 

rf^aV • (41) 

Thus the relaxation rate r-f 1 shows an even stronger quadrupole and magnetic field depen- 
dence than t^ 1 and investigation of dielectric loss or internal friction in the presence of the 
magnetic field opens an attractive opportunity to investigate the role of the R~ 3 interaction. 

The above analysis is based on the assumption that the EFG has approximately axial 
symmetry, i.e., that the parameter x is small. Generally this is not the case. Nevertheless, 
the magnetic field induced relaxation mechanism has a quite universal character. Indeed, 
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let the nuclear spin be half-integer. It follows from Kramers' theorem that in zero magnetic 
field the energy spectrum is degenerate. So, after application of the magnetic field the total 
number of levels of the multilevel TS increases and an effect takes place. 

The nature of the tunneling systems in amorphous solids remains puzzling despite of the 
large theoretical efforts probing various models. The main problem of the theory is the lack 
of experiments which test particular models versus the original phenomenological models, 
which simply employs the distribution (fT|). We expect that the magnetic field experiments 
allow to reveal which and how many atoms participate in the tunneling. This should shed 
new light on the microscopic nature of the tunneling systems in glasses. 
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